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ABSTRACT 


Adaptive arrays are highly versatile sensors employed in modern wireless communication 
systems to combat interference and multi-path fading and thereby increase system capacity. In 
adaptive processing, weights are attached to the incoming signal at each element to produce nulls 
in the directions of interferers. However, mutual coupling is normally ignored in such processing. 
Instead, the principle of pattern multiplication is used, where the assumption is that the radiation 
pattern of an array is the individual radiation pattern of the elements multiplied by an array 
factor. This assumption ignores mutual coupling and the error can be significant under certain 
conditions. This work begins by setting up the scenario of an adaptive array in a mobile 
communications scenario in a scattering environment. Following that, we introduce the theones 
of mutual coupling and beam-forming. Expressions for the optimum solutions and the signal-to- 
interference-plus-noise ratio are then derived based on the preceding discussions. Subsequently, 
we show the effects of mutual coupling and the conditions under which the effect is significant. 
In addition, we will examine various parameters involved in our scenario and how they affect the 


performance of the adaptive array. 
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I. INTRODUCTION 

Array antennas are well known for their versatility in many areas, not least of 
which in point-to-point communication. By a process of beam-forming, an array antenna 
can steer its main beam in the direction of its intended transmitter or receiver while 
simultaneously pointing its nulls to cancel interference sources. Furthermore, adaptive 
algorithms have been introduced [Ref. 1 : p586 - 591] to infer and respond to the 
directions of these (desired and interference) sources, the positions of which change with 
time. There is currently a lot of interest in using adaptive antennas in mobile 
communications to increase systems capacity or reduce interference. 

However, the effect of mutual coupling has received little attention in the adaptive 
array processing community. It is often assumed that the beam pattern of an array is the 
individual pattern of an element multiplied by an array factor, in which mutual coupling 
is ignored. In this work, we will show that this error is significant under certain 
conditions and examine the effect mutual coupling has on the performance of array 
antennas. 

We begin in Chapter II by describing the scenario of mobile communications on 
which this work is based, followed by theoretical discussions on mutual coupling and 
beam-forming. Finally, we shall derive useful expressions for optimum adaptive solutions 
and for the signal-to-interference-plus-noise ratio (SINR), which will serve as the 
yardstick for measuring the performance of an array antenna. 

In Chapter IIL, we will show the effect of mutual coupling and the conditions 
under which this effect 1s significant. In addition, we will examine various parameters 


involved in our scenario and how they affect the performance of the adaptive array. 





II. THEORY 
In this work we will look at the effect of element mutual coupling on the signal- 
to-interference-plus-noise ratio (SINR) of an adaptive array whose design is based on the 
minimum mean square error algorithm. The elements of the array are assumed to be 
dipoles and the analysis will be carried out in terms of the array mutual impedance. A 
variable number of sources of interference as well as variable inter-element spacing will 


be considered in the study. 


À. DESCRIPTION OF THE PROBLEM 

We consider the case of K+1 transmitters (one desired and K interference) all 
radiating narrowband signals to a receiver at the same carrier frequency. The scenario is 
illustrated in Fig 2.1. The transmitters (also referred to as sources) are all embedded in a 
scattering environment. The scatterers, which could represent reflecting obstacles 
surrounding the source, are assumed to be distributed uniformly around each source over 
a circular region of radius ry. We assume single bounce scattering of waves, meaning that 
waves reach the receiver via single scattering off the obstacles. We further assume a 2D 
situation where waves travel only in the azimuthal (horizontal) plane. Because of the 
presence of scatterers, waves will not arrive just from a single direction for each source 
but from an angular region determined by rx and the distance Rx between the k” source 
and the receiver. The receiving antenna is assumed to be an N-element linear array with 
an inter-element spacing d. The elements are assumed to be center-fed half wave dipoles 
oriented along the z axis, so the pattern of an isolated element in the azimuth plane is 


constant. The mean directions of the desired source (k=0) and the k™ interferer (k=1 to K) 


G2 


are denoted by $, and ¢, respectively. Let 0, be the angle from any scatterer around the k" 


source to the array and 0% = ( - ¢,, as shown. The angle of arrival q is then measured 


from the mean direction @, for each source. 
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Fig 2.1 Wireless Communication Scenano 


Because of the multi-path phenomenon, the amplitude, phase and angle of arrival 
of waves from each source will be random variables. The interferers are assumed to be 
arbitrarily located about the receiver. The purpose of the array 1s two-fold: 

1) the array weights are chosen to produce nulls in the direction of the 

interferers, and thereby reduce their influence on the received signal, and 


11) if N > K+1, to provide diversity gain against multi-path fading. 


The weights of the array are continually changed depending on the location of the 
interferers. Such an array is termed as adaptive array or smart antenna. The situation 
considered above may well pertain to the case of wireless network in a tactical battlefield 
or to the case of mobile cellular radio where K+1 co-channel mobiles send signals to a 
base-station receiver. We have assumed thus far that the mean directions of the desired 
source and the interferers are known a priori. This is seldom the case in reality. In a 
practical mobile situation, the sources move around in space and their directions are not 
known precisely. In such a case, the weights are determined by means of either the Least 
Mean Square (LMS) or Recursive Least Square (RLS) algorithm. However the present 
case is still of interest in applications such as wireless local loops where wireless service 
is provided to fixed subscribers and the directions of subscribers available from a 
previously compiled data base. 

Since the calculation of optimum weights requires the covariance matrix of the 
received signal vector, we would like to know the probability density function (pdf) of the 
angle of arrival. For a uniform density of scatterers about the source, the pdf of the angle 


of arrival with respect to the mean direction is given as [Ref. 2 : p24-25] : 


2 2.2 
p(O&) = —. cosa, V1 - bj sin' a la, |< m (1) 


where Dk = Rk/Tk (2) 
and C, is the maximum angle of arrival from the k™ scattering region, given by 
bisinCk = | 
For subsequent analysis, it is simpler to express the pdf in terms of a new angle Yk 


defined by siny, = bysino ; then 


cos Yk dk = by cos oy doy, (3) 


Using this transformation of variable [Ref. 3 : p122], the pdf of y 1s: 
da, 
dy x 


2 
= cos” yy - sp <= (4) 
T 2 


pO) == P(Ox) 





This will be applied in Section 1D when the expressions for optimum weights and SINR 


are derived. 
The mean value of O& 1s @,= 0 and its variance à for by»»1 1s approximately 


given as [Ref. 4] : 


= 3 
a =b- (5) 


2 
T A (| 
2 
2b, 


The root mean square (rms) value of the angular spread is then ok For example, with ry 


ü 


= 100m ana, — | knee, = 5.7 and Ot. 285. 


B. MUTUAL COUPLING 

In an array of electrically identical antenna elements, it is common to assume the 
total radiation pattern to be the product of the array factor and the radiation pattern of the 
individual elements, i.e., total radiation pattern = (array factor) x (element radiation 


pattern). Implicit in this formulation is the assumption that the individual radiation 





patterns remain the same even when placed in proximity of each other when mutual 
electromagnetic interaction exists. However, such an assumption is erroneous in view of 
element mutual coupling. When an antenna array is designed to provide diversity gain, 
mutual coupling will disturb the spatial correlation of signals received by the array. In this 
thesis, we will examine the effect of mutual coupling on the performance of the array. 

The input impedance of an isolated antenna is simply the ratio of the voltage 
applied (at its terminals) to the resulting terminal current. This is known as the self- 
impedance and is, in general, a complex value as both voltage and current are of different 
phases. When two antennas are located in close proximity, the current distribution (and 
hence the radiation pattern) on one is affected by the field radiated from the other. This is 
known as mutual coupling. To account for the mutual coupling effect, extra impedance 
terms called mutual impedances are needed. 

In general, the solution to the mutual impedance for dipoles can be accomplished 
by way of coupled integral equations. Apart from its dependence on physical dimensions 
and separation distance between the elements, the mutual impedance also changes with 
the number of antenna elements. For an array of N elements, the impedance matrix can be 
written as Zen denoting the mutual impedance on the m? element induced by the pn? 
element in an N-element array. The details of mutual impedance calculation are not 
within the scope of this paper; we are interested in the effects it has on array antennas. We 
will use results from simulation programs or from analytical expressions. As mentioned 
earlier, for the sake of simplicity, we shall only concern ourselves with half wavelength 
dipoles of small thickness. Even though Zmn .! varies with N, we will assume Z4, ! — 


Ze actual computation of Zu and Zaun show them to be very close to Zo (The 


superscript I is subsequently dropped.) Zmn, the mutual impedance between two side-by- 
side dipoles, each of length L and separated by a distance d, is given by [Ref. 5: p265] : 
Zi; - 30(2E; (-jkod) - E;[- jk, (Vd? € L2? - L] - E;j[- ko (vd? - L? - L)]) (6) 

where k, 2 21/X is the wavenumber in free space, and 

E(t jy) = Ci(y) € j Si(y) 

Ci(y) 2 cosine integral 

Si(y) = sine integral 

It is obvious that the inter-element spacing d has an influence on Z4; as the 
distance apart increases, one can expect the contribution from other array elements to be 
reduced and hence Zp, (m #n) diminishes. The resultant impedance matrix Z is then 


almost diagonal and mutual coupling can be ignored. 


E. BEAM-FORMING ALGORITHM 
Consider an arrangement of N elements spaced d apart (Fig 2.2). An 
electromagnetic wave arriving from bearing 6 induces an open circuited voltage V,” on 


" element. The open circuited voltage depends on the incident field and the 


the n' 
orientation of the dipoles with respect to the incoming wave. The coupling between the 
elements is governed by the impedance matrix Z. It 1s assumed that the adaptive 
processor presents load impedance (Zi), (nz1, 2, ... N) at each element. The terminal 
voltage and current at the n port are denoted by V, and I, respectively. It is clear that V= 


-Zil, where V is a vector of voltages presented to the adaptive processor and Zy is the 


diagonal matrix 2 diag(Zi1 , Zi» , ... , Zin). 
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Adaptive Processor - 


Fig 2.2 Array of N elements 


The terminal voltage is related to the open circuited voltage by means of the 
impedance matrix as : 
Vi = Zul + Zib + Zib +... + Zınin + Vi 
Və = Zal + Zao + Zy3]5 + ... + ZonIn + V2 


V3 = Zaılı sp Zyb Al 23313 gum ata Z3NIN + VO 


Vn = Znılı + Zub + ZN313 ee eee ats ZNNIN + Vn 
or equivalently, 


V =ZI + V (7) 


For a dipole antenna the open circuited voltage can be obtained from the effective 


length of the antenna by 


Eos = h a ER (8) 
where h = vector effective length of n” element 
EYS = incident electric field on n element. 


For dipoles of length L, oriented along the z-axis and having a sinusoidal 
distribution I(z) = I, sin{k,(L/2-Izl)}, the effective length is 


f(0) 




















xs (9) 
T en 
sın 
2 
where 0 1s the usual polar coordinate in the spherical coordinate system and 
k IE 
cosi > cos 8) — cos( > ) 
f(0) = E 
> sin O 
is the radiation pattern of the dipole. 
The electric field incident on the n^ dipole is of the form, 
Eius B GB) E, 
In the present case, 0 — 7/2 and the radiation pattern of the dipole reduces to 
f(0) = (1 - cos Kol y 
2 
ke 
,  l-cos( 9—) | | 
Therefore, V? = ČE, uec ae = Ve jsdin-Dase (10) 
: O 
sin 
( E ) 
E 
where V, 2 —P-tan(——) 
4 
Using (10) in (7) and substituting V = -Z, I , we get 
V UNS eu 





where LA LS (12) 
and Zo = (Z + Z) (13) 

It is this voltage vector V and not V that is available for array processing. It is 
seen that the effect of mutual coupling is to couple the open circuited voltages induced at 
each element. Thus mutual coupling will affect the spatial correlation between the 
element signals and hence affect the diversity gain of the antenna. 


In the presence of thermal noise, the total voltage input to the adaptive processor 


V=ZV”+N (14) 
where N(t) = column vector of thermal noise. 
With this input signal vector V, , the adaptive processor attaches a weight w, (* is 
the complex conjugate operator) to each (Vi), as shown in Fig 2.3, and produces a 


combined output S,. 





Fig 2.3 Adaptive Processor 


In mathematical terms, the array response is 
* * * 
S;=W, Vi +w2 V>3+....+Wnu VN 


=w" V, (15) 


where " denotes the Hermitian operator (transpose and conjugate). 
By varying the weighting vector w, up to N-1 nulls can be steered in any direction. 
For example, in the absence of mutual coupling and scattering, and using a 4-element 


array, if the nulls are desired at 4, , d2 and $5 we choose the weights according to 


M: 
V v -] 
VS m (16) 
Ve ae m 
3 3 3 W, 


Wi | 
where Ve = e ort eos Oy 1,93. Ford) = 45% 05 = 90°, 0, = 135% nd d ee 


we get 2 =0.2114, 23 =-0.2114 and =1. The array response to a plane wave 
Wi W, Wi 


arriving from angle 6 is plotted in Fig 2.4. 
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Fig 2.4 Beam pattern of a 4 element array, intended nulls at d = 45°, 90° and 135? 
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IU: Criteria for Optimum Weight 

It was shown in the preceding paragraphs that, by using an appropriate weighting 
vector, one is able to direct the nulls towards the direction of interference sources, therebv 
cancelling them out. How does one choose the weights in the presence of noise and 
scatterers? 

The choice of the weighting function is based on the statistics of the total signal 
received (desired + interference). Basically, the objective 1s to optimize the beam towards 
the desired signal according to a prescribed criterion, of which there are several. It has 
been shown [Ref. 5: p41] that the various criteria (including the Minimum Mean Square 
Error (MMSE) and the Maximum SINR criteria) converge to the same Wiener solution, 
differing only by a trivial scalar factor. Since the solutions are all the same, we shall use 
the MMSE cniterion. In the MMSE algorithm, the weights are chosen to minimize the 
error, defined as: 

error = <| A,(t) — S,(t) |? > (17) 
where < > is the expectation operator. The weights that minimize the error are shown to 


be of the form [Ref. 6: p39] : 


Wopt = ® S (18) 
where & = covariance matrix =<V, Vi (19) 
S = <A," V> (20) 


A, = reference signal = ,exp[j(@t + Bo)] 
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D. Choice of inter-element spacing d 
Before proceeding to the calculation of the weight vectors, it would be instructive 
to consider the effect the separation distance between array elements on the input signal. 


Suppose there were no mutual coupling, no thermal noise (N(t)=0) and all loads are 


, | Z 
identical (hence Z. = z.U , where U is an NxN identity matrix and z, = m Then 
L t41] 


using (10) - (14), the input signal to the array processor is : 


(Von N Es A (21) 


Suppose the weights are all equal, i.e. W, = w for all n, then 


S(0 - Y w,*V, 


N 
jkod(n-1)cos 9 
WI O y eo 
l 


_ a JNK dcos > 
l e (22) 


w *ZeVo —————--— 
0 1 eJKodcos à 


The maxima are obtained at cos 6 = + n(4), (n=0,1,2...) where the denominator 
d 


becomes zero. In general, there is the desired major lobe and other major lobes called 


grating lobes. To avoid these grating lobes, one chooses 4 > 1, so that the above equation 
d 


will have only one maximum corresponding to n=0. Therefore in this paper, we shall be 


primarily considering values of d up to IA for most results. 
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The nulls of S,(@) occur when the numerator alone goes to zero, 1.e., whenever 


gJNK,dcos 9 = 1, assuming k d cosd #2mn. This means cos 6 = + ee , where m is 
Nd 
any integer. As such, for 4 <1, there are several solutions for ò and m, which implies the 
d 


presence of multiple nulls, which might dramatically affect the SINR if they coincide 
with the direction of the desired source. The above example assumed equal weights and 
the absence of noise and mutual coupling. If these were taken into consideration, the extra 


lobes would, in general still be present, albeit shifted. 


D. EXPRESSIONS FOR W,,: AND SINR 

With the results of the first three sections, we are now ready to proceed to derive 
the expressions for the optimum weighting vector (Wopt) and the signal-to-interference- 
plus-noise ratio (SINR). 

To begin with, let the desired signal be of the form A,(t) V4** and the interference 


signal be of the form A,(t) Vx“ where 


A(t) = a exp[j(at + Bo)] (23) 
Ax(t) = ax exp[j(a@t + BoJ (24) 
(Vor In = Vo explj(ko d (n-1) cos 60)] (25) 
(Vi )n = Vo explj(ko d (n-1) cos $)] (26) 


Ay , Ax = amplitudes of signals 


Bo , Bx = phases of signals 
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Then the input signal to the adaptive processor is 


V,=Z.V” +Ntt) 


K 
= Z AA Vo + Y Ax Vio) +N(0) (27) 


k=1 
I Optimum Weight Wop 
We now attempt to find a convenient expression for the covariance matrix ®. 
Using equation (19) : 
pz2«V, VB» 
-«(Z. V* e N) (Z. V* - N)! » 
22.« Velas NS 
+ Z.<V“NHS+<NN"> 
The following assumptions are made : 
1) The noise voltage at each port ıs independent of the each other and of the 
signals; its pdf is assumed to be zero mean Gaussian. 
ii) The open circuited desired signal is assumed to be independent of the open 
circuited interfering signals. 
111) The open circuited interfering signals are assumed to be independent of 
one another. 
1v) The amplitude, phase and angle of arrival of the signal from each source 
vary independently of each other. 
It is shown in Appendix A that this eventually reduces to : 


DEV UE C LIA OZ Ex Cua Los 


VA UE ICAA DA LACAN (28) 





where Va = mean square noise voltage at each port 


U is the N x N identity matrix 


Z.— as deimedan eqn (2) 





pM ODE EI T) 


kgd(m - n)sin OTI 


Ry 
| | ME 
čo = input signal (desired source) to thermal noise ratio = ——————— 
MS 
th: Ne < 2 > 
Ek = input interference (k interfering source) to thermal noise ratio = -p 
VN 
J; (.) = Bessel’s function of the first kind of order one 
It is also shown in Appendix A that : 
S= <A Vo 
= Vo es Le e (30) 
where 
2]; k.d(m - sin dur, 
[Cy m= eJkod(m—l) cos Do Ro = (31) 
k,d(m — 1) sin ofo 
Ro 
Substituting eqns (30) and (31) into eqn (18) : 
Wopt = o! S 
= ay Eo [U + ZeGoCo Ze + Ze E EkCr Ze” J” Z. Cs (32) 


where a, € 1s a trivial constant factor and is subsequently dropped. 


2 Signal-to-Interference-Plus-Noise Ratio 
With the weight thus calculated, we can proceed to evaluate the SINR, which 


gives a measure of the effectiveness of the array : 


SINR = Power Received from desired signal 
Noise Power + Interference Power 


Po 


a ae (33) 
Py + >) Py 
k=) 
where, from Appendix B, 
Po=<(w' Ze Ag Vor) (W' Ze Ao Von)" > 
=V? Ew" Ze Co Ze" W (34) 
and 
P< (wW Z A VEIN Z Ar Veo 
= V? w" Ze Ek Ck Ze” w (35) 
Pu = < (w" N) (w” NË >= VI ww (36) 
H H 
Z Z 
Therefore, SIN A olan (37) 


wHw+wH Lame er ZE y 


III. RESULTS, OBSERVATIONS AND DISCUSSIONS 
This chapter discusses the results obtained and presents explanations of these 
observations. We first start by asking whether mutual coupling is at all significant and 
then proceed to look at how separation distance, noise level, scatterer radius and number 


of interference sources affect the performance of the array. 


A. SPATIAL CROSS CORRELATION 
The first question we need to ask is whether the effect of mutual coupling is at all 
significant to an adaptive array? The importance of mutual coupling to diversity arrays 
has been addressed in [Ref. 7]. The parameter of interest is the spatial correlation of the 
field. It is an indication of the correlation of the voltage received between the elements of 
an array. One of the benefits of an array is to provide diversity; this means that while one 
element may be in the deep fade of the received signal pattern, hopefully another would 
still be in a position to receive the signal. For this to be true, the elements have to be 
relatively uncorrelated — hence a small value of spatial correlation means greater 
diversity. 
Consider a two element array, distance d apart. In the absence of thermal noise, 
the voltage equations from (14) becomes : 
Vi = Zen Vi” + Ze n V2” 
V2 = Zeca Vi" + Zen V2” 
Then the spatial cross correlation of the received electric field is defined as : 
p =< Vi V? >= Zen Zer <V” V” > + Zen Zez <Vi% Vo" > 


* * * * 
+ Zc Zez <V2 V >+ Zeiz Zez «V^ V > (38) 





MY = ao c (39) 
and «VIT Vo" > =<¥> Vi >= Mi (Ref. 4] (40) 
where Zen = Zez , Zen = cz (41) 

Jo (.) = Bessel’s function of first kind, order zero. 
Thus 

p = (Zen Zeiz? + Zen Zen )+ (Zen Zen + Ze Zen) (kd) (42) 

Extracting the common factor and substituting for Zc, the spatial cross correlation 
of the electric field for a two element array terminated in Zi, = Zy is given by [Ref. 4] : 
R, 


P = Po - — (43) 
X ler 
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where Pæ = Jo(kod) is the correlation in the absence of mutual coupling, and Ry; and Rj» 
are the real parts of Zi; and Zi respectively. Fig 3.1 below shows the spatial cross 


correlation vs separation distance d between two half-wavelength dipoles. 


Spatial Cross Correlation of a 2-element array 
0.35 


— — With Mutual Couplin 
——— Without Mutual Coupling 


0.3 


0.25 
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I< VV; >l 





distance apart d/A 


Fig 3.1 Spatial Cross Correlation of a 2 element array 
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As expected, at longer separation, the effect of mutual coupling is diminished and 
the two curves converge. However, there is a marked difference for d<A, which indicates 
the significance of including mutual coupling. 

Another way to interpret spatial cross correlation is to examine the normalized 
spatial cross correlation (Fig 3.2). For reference, we have indicated the correlation values 
of 0.7 and 0.3 (Table 3.1). It is seen that for this example the effect of mutual coupling is 


to decorrelate the signals initially but correlate them more for larger values of spacing d. 
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Fig 3.2 Normalized Spatial Cross Correlation Electric Field in a 2 element array 







Values of d/A 


Without mutual coupling With mutual coupling 
2 1.1201 2 3.0701 
p € 0.7 2 0.1801 2 0.0901 


Table 3.1 Values of spacing d for p € 0.3 and 0.7 
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While both curves fall rapidly below the upper level of 0.7, there is a large 
difference in the lower level of 0.3; without mutual impedance, the arrays have minimal 
spatial cross correlation (<0.3) for d>1.12A, whereas with mutual coupling, this occurs 
only after d>3.07A. What this means is that in the presence of mutual coupling, the array 


diversity will be generally reduced. 


B. EFFECT OF MUTUAL COUPLING ON WEIGHTING VECTOR 
Another effect of mutual coupling is in the resultant weights. Using (18), we 
Observe that we get different weighting functions when considering and ignoring the 


effects of mutual impedance (Table 3.2). The weights are normalized such that w” w = 1. 





Number of elements N=4 
Separation between elements d=)/4 
Number of interference sources K=3 
Direction of desired source ¢, =90° 
Directions of interference sources é, _ 30°. 60°. 120° 


Scatter radü Or 
Consider mutual coupling Deterministic case 
w (normalized) w,(normalized) | wa(normalized) 

0.2624-0.1203i 






















0.4330+0.56001 0.2335+0.54451 0.4275+0.45641 
0.1704-0.35301 -0.0689-0.32271 


0.2102-0.32111 


Table 3.2 Weights for adaptive MMSE and deterministic cases 


Table 3.2 also shows the weights wa calculated in the deterministic case, that is, 
where it is assumed that the directions of the interference sources are known and the 


weights calculated such that nulls are steered in the directions of the interferences. Since 
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the weighting vectors are all different, we would expect different beam patterns using 
these 3 solutions, as shown in Fig 3.3 below. Values of the corresponding SINR are 


indicated in the inset. 


Note : In Fig 3.3 and subsequent plots, the blue curve corresponds to adaptive (MMSE) 
solution using eqn (18) , with mutual coupling considered. This means that the off 
diagonal terms of Z are non-zero and are given by (6). The red curve corresponds to 
adaptive solution also given by (18), this time ignoring mutual coupling (Z in (13) is 
diagonal). The green curve corresponds to the solution in the deterministic case, 
calculated without mutual coupling. Moreover, the direction of the desired source is 


maintained at @ =90° (unless otherwise specified) and is denoted by a blank circle on the 


abcissa. Pink asterisks (*) denote the directions of the interference sources. 


Weights normalised. N= 4, d= 0.254, K =3, r /R=0, 1, /R=0, E, = 20dB, E, = 20dB 
1 


0.9 —— With MC, SINR = 8.6dB 
“~~ Ignonng MC, SINRo = 6.148 


terministc, SINRd = 3.0dB 
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Fig 3.3 Beam pattern for adaptive solutions considering and 
ignoring mutual impedance and deterministic weights. 
All weights normalized such that w'' w=1, beam pattern not normalized. 
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In Fig 3.3, it is seen that although nulls are steered in the directions of the 
interference sources in the deterministic case (green), it does not provide the best SINR. 
This is because in the other cases (adaptive solutions — blue and red), the gain in the 
direction of the desired signal is higher than the deterministic case (which does not take 
into account direction of desired signal), thereby contributing to a higher desired signal 
power received, P, in (34). This thus gives rise to higher SINR. 

Fig 3.3 shows the beam pattern (un-normalized) using weights calculated. 
However, it is often the normalized beam pattern that is more suitable for comparison 
between the three solutions; in Fig 3.4 and subsequent graphs, the beam patterns are 


normalized to the maximum gain. 
N= 4, d= 0.254, K = 3, 1, /R=0,1,/R =0,& =20dB, &, = 20dB 


09 ——— . With MC, SINR - 8.6dB 
: — Ignoring MC, SINRo = 6.1dB 
ee eterministc, SINRd = 3.0dB 
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Fig 3.4 Beam pattern (normalized to max gain) tor adaptive solutions considering and 
ignoring mutual impedance and using deterministic weights 


Fig 3.3 and Fig 3.4 are essentially the same beam patterns, though in Fig 3.4 the 


adaptive, no-mutual-coupling solution (red) seems to have a higher gain in the direction of 


the desired source (¢, = 90°) than the mutual coupling case (blue). As such, one might be 
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lulled into expecting the red curve to have a higher SINR; it should be borne in mind that 
the graph shows relative gains only. For a true value of the SINR values, one should 
always refer to the values displayed in the inset rather than to the heights of the curves. 
The previous two sections show that mutual coupling does have an effect on the 
performance of adaptive array antennas. This effect obviously diminishes as the physical 


separation between the clements increases. 


C. EFFECT OF ELEMENT SEPARATION 
An interesting comparison would be how the SINR is affected by varying the 
separation distance d between elements. Fig 3.5 below shows the SINR vs d plot for a 4- 


element array antenna in a 3 interference situation (no scatterers). 


N= 4, K=3, 1, /R=0,1,/R=0,£,=20dB, £, =20dB. q, = (30%, 60”, 120”) 


SINR (dB) 
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Fig 3.5 SINR vsd 
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It is seen that the SINRs with and without mutual coupling do show differences for 
small values of d (less than 1A), sometimes by as much as 3dB (e.g. at d=0.75A). At 
higher separations, they are about the same, as expected since mutual coupling effect is 
diminished. In all cases, the SINR using the adaptive solution (no mutual coupling) is 
consistently higher than that using the deterministic weights. At the same time, there are 
differences between the adaptive solutions with and without mutual coupling (blue and red 
curves in Fig 3.5). The exact difference is difficult to quantify as it is case specific: in Fig 
3.5 above, the situation is for a 4 clement array in a 3 interference environment (Ó, = 30. 
60° and 120°) with & = & = 20dB. 

An interesting feature is the fact that the SINR drops drastically near d=1.154. 


This is shown in Fig 3.6 below. 
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Fig 3.6 Beam pattern for d = 1.15A 


A closer examination reveals that this drop is due to the rather low angular 
response of the array in the direction of the main source. Therefore, in practical design of 


array antennas, it would be wise not to exceed d=A/2 unless we impose an additional 
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constraint that the pattern not go to zero in the direction of the desired source. In this 


paper, we shall continue to study the cases for d up to IA. 


D. EFFECT OF NOISE 
From (18) and (14), we have : 
Wop = 9” S 
-[U+Z EG ZIIT 2 E 7 C 
and V, 2 Z. V* + Nit) 


In the situation where & >> 1, & >> 1, we get IV**l 2» IN(t)I so 


Wopi = [Ze (Eo Co) Zo” + Ze (EEx Ca) Ze Y” Zo Cs (44) 
and Mis Lave (45) 
The resultant signal 1s 

Sr = Wopt Vi 


-([Z&C.Z' -Zr&5cz']zcyzv" 
202 7. 17.5, Cu 220, 7 21002 

C6 ZU &SCHZ LE ZN ZN 
=C ZE (Z) (6 Got De Cou ale Zs ve 
=C." (Ze"Ze)" (Go Co + Z & Cy [Ze rv 


= (ey Cae > CC elem (46) 
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which turns out to be independent of Ze. Since the effect of mutual interference is in the 
impedance matrix Ze, it is predicted that when &, , & is large (240dB), the effect of mutual 
impedance is reduced and the two adaptive solutions (blue and red in Fig 3.7) converge. 

(Recall that &,= input signal (desired source) to thermal noise ratio and & = input 


interference (k" interfering source) to thermal noise ratio ) 
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Fig 3.7 Comparison of the effect of noise 


Fig 3.7a shows the beam patterns using the different weights when €, = & = 20dB. 
However, when signal-to-noise-ratio &, and interference-to-noise-ratio €; are increased (to 
40dB) as in Fig 3.7b, one can casily observe that the solutions (with and without mutual 
coupling) converge, as predicted. It is important to note that though the normalized beam 
patterns in Fig 3.7b are similar, the weights and terminal voltages with and without mutual 
coupling are not the same. This is evident in the fact that the SINR values with and 
without mutual coupling differ by as much as 4dB in Fig 3.7b. This is because the noise 


levels in the three cases are different relative to the output signal. 
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E. EFFECT OF SCATTERER RADIL r, , ry 

Recall that in our model outlined in Section 2A that the desired and interference 
sources are at distance Ry from the array. In addition, each source is embedded in a 
scattering environment, in which the scatterers are assumed to be distributed uniformly 


around each source over a circular region of radius r, (k=0,1,....K). Here we take Ry = 1 
km and vary the values of rm (k=0,1,....K); ultimately, it is the ratio b, = that is 


significant in the calculations. 

It was found that the effect of the scatterer radii ro and rg depend on the input 
signal to noise ratios € and & . First we compare the effect of increasing scatterer radius 
To While maintaining ry = 0. For illustration, a 4-element array with A/2 spacing is used in a 
3 interference situation (d, = 30°, 45°, 60°). Fig 3.8 below shows the resultant normalized 


beam pattern. 
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In Fig 3.8a, there are no scatterers around the desired source (r, = 0) whereas in 
Fig3.8b, scatterers are present (r/R = 0.1). We can see that the weights and (hence) beam 
patterns are changed slightly, and that with scatterers (Fig 3.8b) the resultant SINR values 
are slightly reduced. The light blue column shows the extent of the angular spread of the 
desired signals; hence the area in the column under each curve would represent the 
amount of power arriving from the source. In the calculation of SINR, this area is divided 
by the angular spread to give an average value of the signal. Whether this average value is 
greater or smaller than the SINR in the no-scatterer case obviously depends on the shape 
of the curve. In Fig 3.8 above, it is observed that though the average SINR values for 
adaptive solutions are smaller, the deterministic solution actually shows an improvement in 
the averaged case. 

Nonetheless, it is important to note that the adaptive solutions (blue and red) still 
give greater SINR than the deterministic case (green). This is in accordance with our 
conjecture that the adaptive MMSE solution also leads to the Maximum SINR solution. It 
has been shown in [Ref. 4] that the SINR without mutual coupling and no scatterers 
around the desired source is always better than the deterministic case. No such conclusion 
can be drawn when mutual coupling is taken into account. 

Next we investigate the situation when €, and & are increased to 20dB (Fig 3.9); 
comparing the 4 cases (in Fig 3.8 and Fig 3.9), the deterministic case remains the same 
since the deterministic weights do not depend on r, , r; , & nor ék. The adaptive solutions 
change accordingly - most surprising is the fact that for r/R — 0.1, n/R 2-0, $;, 2 &k = 


20dB (Fig 3.9b), the adaptive solutions (both with and without mutual coupling) give 
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SINR values smaller than the deterministic case. This is in contrast to the no-scatterer case 


where the adaptive (MMSE) solutions always give the maximum SINR. 


N= 4, d=0.51,K=3,r /R=0,r,/R=0,E =20dB, £, = 20dB 
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Fig 3.9 Effect of scatterer radius r, (€ and & = 20dB) 


Next we examine the effect of increasing the scatterer radius rą while maintaining 
ro; =0. In Fig 3.10 below, the columns in magenta represent the angular spread of the 
interference signal from the interference sources. The increment of rm again reduces the 
SINR but not as much as in the case of scatterers around the desired source. The decrease 
in SINR is expected since scatterers increase the directions from which the interference 
signals can now reach the array, thereby increasing the interference power received. In this 
case the increment of ry does not change the fact that adaptive solutions still give greater 


SINR than the deterministic case. 
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Fig 3.10 Effect of scatterer radius r, (& and & = 20dB) 


F. EFFECT OF NUMBER OF INTERFERENCE SOURCES K 
Another variable parameter is K, the number of interference sources. In Fig 3.11 
below, a 4-element array (N=4) is placed in an environment where there are no scatterers 


around any of the sources (ro = rk = 0). 
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Fig 3.11 Effect of Number of Interference Sources K (K>N-1) 
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In Fig 3.11a, 3 interference sources (K=3) have been specified in the directions 
Ó, = 30°, 60°, 120°. The adaptive solutions and the deterministic solution are similar and all 


steer nulls in the directions 30°, 60°, 120° that effectively cancel out the interference 
sources. The resultant SINRs for all 3 cases are therefore very close. In Fig 3.11b, an 


additional interference source has been added (at ¢,= 140°). With the deterministic 


weights, a maximum of N-1 nulls can be chosen; this means we can (arbitrarily) choose 
only 3 of the 4 interference sources to cancel. With the adaptive solutions specified by 
(18), the weights are altered to give a different beam pattern as shown. Though none of 
the nulls is exactly in the direction of the interference sources, the ‘spreading’ of the 
pattern is such that the resultant SINR is greater than that given by the deterministic case. 
In Fig3.12 below, the same situations are compared, this time with scatterers 


around all sources. 
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Fig 3.12 Effect of Number of Interference Sources K (K>N-1) 
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In Fig3.12a, we note that the deterministic SINR is better than for the adaptive 


cases (as discussed previously in Section E). However, with the addition of a fourth 
interferer at db, = 140°, the adaptive cases actually perform better than the deterministic 


solution in terms of SINR. This emphasizes the flexibility of the adaptive solution in 
encountering more than N-1 inteference sources. 

At the other extreme, we would like to examine the situations where the number of 
inteferers K is less than the number of array elements N. To illustrate, we consider N=6 
and 8 in an environment with only 3 interference sources. Fig 3.13 below shows the 


resultant normalized beam patterns. 
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Fig 3.13 Effect of Number of Interference Sources K (K<N-1) 


It is seen that when K « N-1 ("there are more than enough array elements to go 
around"), the SINR using the adaptive solutions are improved. Table 3.3 below shows 


these values for N=4,6,8,10 and 30 (K=3 in all cases). 
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Table 3.3 Comparison of SINR 





SINR (with mutual coupling) 












SINR (ignore mutual coupling) 
SINR (deterministic case based on N=4) 













As is obvious from the table, as the number of array elements (N) increases, the 
SINR increases as well for the adaptive cases. Of course, the deterministic case does not 
change since it is dependent only on the direction of the interferers. Initially, we note that 
SINR (deterministic) is greater than SINR(adaptive) — this was discussed in Section 3E as 
being due to scatterers around the desired source. However, for higher values of N (N2 8 
in this example), we find that the SINR of the adaptive solutions actually increases to 
greater than that offered by the deterministic case. Beyond a certain value of N, this 
increment becomes insignificant or non-efficient (c.f. SINR for N=10 and 30). 
Nonetheless, this is another indication of the ability of the adaptive algorithm in making 


full use of its elements to find the best solution. 


G. ANGULAR RESOLUTION OF SOURCES 
Our next area of interest is in examining if the use of an adaptive array is effective 
when the various sources are closely spaced. Consider the situation where one of the 


interferers is very close to the desired source : 6, = 90°, ¢, = 90°-2.85° = 87.15° (2.85° 
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ean po : 
was chosen as it is the rms value of angular spread for E =— =().1 as in eqn (5) ). This 





is shown in Fig 3.14. 
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3.14 Effect of Interference Source close to Desired Source (& and & = 10dB) 


Fig 3.14a shows the situation where there are no scatterers around any of the 
sources. It is seen that the close proximity of a null near the desired source lowers the 
SINR in the deterministic case drastically, whereas for the adaptive solutions, the nulls are 
actually steered away so as to give a significantly better SINR. In Fig 3.14b, the presence 
of scatterers does affect the performance of the adaptive solutions, but it is seen that the 
adaptive solutions still give better SINR values. 

Fig3.15a and b below shows the same situations with & and & increased to 40dB. 
As discussed in Section D, high values of €, and € mean that the solutions all converge to 
the deterministic case. However, in the presence of interference close to the desired 


source, the solutions are different. More significantly, the SINR from the adaptive 
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solutions in this case is lower than that from the deterministic case — again, this is due to 
the presence of scatterers around the desired source. 
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3.15 Effect of Interference Source close to Desired Source (&, and & = 40dB) 


What happens when an interference source, initially remote from the desired 


source, is moved close to it? In Fig 3.16a below, $, = 30°, 45° and 60° while in Fig 3.16b, 


the interference source at Ø, = 45° is moved to 87.15”. 
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Fig 3.16 Effect of Interference Source close to Desired Source (&, and & = 20dB) 


In a previous section, it was shown that when there are scatterers around the 
desired source, the deterministic solution gives a greater SINR than the adaptive solution, 
as shown in Fig 3.16a. However, this is reversed when one of the interference sources is 
moved close to the desired (Fig 3.16b). This shows that when an interference source is 
close to the desired source, the adaptive solution performs better than the deterministic 


one, even when there are scatterers around the desired source. 
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IV SUMMARY AND CONCLUSIONS 

A. SUMMARY 

Through the sertes of results obtained, the effect of mutual coupling on the 
performance of adaptive array antennas was studied by varying several parameters. These 
can be summarized as follows: 

1. Effect of Mutual Coupling on Spatial Cross Correlation 

The effect of mutual coupling is to reduce the spatial cross correlation initially 
(for d < 1A) but to increase it for larger values of inter-element spacing d. 

2: Effect of Separation Distance d 

For d»1A, there is practically no difference on the SINR obtained with and 
without considering mutual coupling. At lesser separation, however, the difference is 
significant and can be as much as 3dB. 

3 Effect of Noise 

When the signal-to-noise-ratios €) and & are high (>40dB), it was found that the 
solutions obtained with and without considering mutual coupling are very close, and 
converge to the deterministic solution (provided number of interference K < number of 
array elements N). 

4. Effect of Scatterer Radii r, and r; 

When the scatterer radii ro and ry are increased, in general, the SINR of the 
solutions are reduced, This is expected as it represents uncertainty in the directions of 


arrival of the sources. 
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When the scatterer radius ro is increased, at a critical &, & level, the adaptive 
MMSE solutions cease to give Maximum SINR. This is in contradiction to the conjecture 
that MMSE and Maximum SINR solutions are equivalent [Ref. 6: p41]. 

On the other hand, increasing the scatterer radius ry; alone does not ‘invert’ this 
trend of MMSE solution giving max SINR. 

3. Effect of Number of Interference Sources K 

When the number of interferers K 2 the number of array elements N, the 
deterministic solution can do no better than pointing N-1 nulls towards the interferers. 
The adaptive solutions, however, seeks to minimize the error and generally gives a 
greater SINR than the deterministic case. 

When the number of interferers K < the number of array elements N, it was found 
that the adaptive algorithm makes use of this fact to improve the SINR. Hence even with 
scatterers around the desired source (as in the previous section), with a sufficient number 
of array elements, the adaptive solution will give an SINR greater than that in the 
deterministic case. 

6. Angular Resolution of Sources 

In cases where one of the interferers is very close to the desired source, it was 
found that adaptive solutions perform better than deterministic one. This includes the 


situation where there are scatterers around the desired source and demonstrates the 


excellent angular resolution offered by adaptive arrays. 
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B. CONCLUSIONS 

We have examined the effect of mutual coupling on the performance of an 
adaptive array in a realistic communications environment, where the desired and 
interference sources are embedded in an area of scattering objects. Expressions for the 
MMSE solutions and SINR were derived for this scenario. These were subsequently 
used to show the effect of mutual coupling, in particular in array where the element 
Separation was less than 1A. The effects of the other parameters in the system were also 
discussed, including conditions under which the MMSE solution does not give rise to 


Maximum SINR. 
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APPENDIX A : EXPRESSIONS FOR ® AND S 


A. EXPRESSION FOR ® 


then 


SO 


From (15), 
covariance matrix © = <V, V> 

= < (Z. V“ + N) (Z. V + NË > 

=Ze< VV E> Zo + Ze< VON > (Al) 

EEN VES E NINA > 
Now, consider the terms separately : 
(1) Second and third terms : 
Assuming noise is zero mean and uncorrelated with signal, 
ds e XQ» B) Nt) > ZN) >=0 
(A2) 
< VEEN” >=<V><N"> =20=<NVrH> (A3) 
(11) Fourth term : 
Assuming noises at the ports are uncorrelated with each other 
< N NË >mn = ôm VA 


(A4) 


where Vy = root mean square noise voltage 


(11) — First term: 
= Vo 2 H an [as e)(wt+Bo) vs qd eg (GU Bo Wace ] 
| [ag e (QU Bo Mas H + y ay e XQ bb NDS H ] > 


Sian Vio (WV) Se = aya ORS A a 


43 


a es OVE ee el ee Wace (Vi eS 


Now a, Vi M c Mu (A6) 
(assuming magnitudes of narrowband signals are uncorrelated to the angles of arrivals) 
Also, assuming the interfering signals are uncorrelated with the desired signal, 
then. <a oak = <> 0 
since the phase Bo is assumed uniform over (0, 2x). So 
<a ay e PEPE NAS = > ca en oC E 
=0 
and. <a ni 


ey > NER (A8) 


Now consider Ni (a e e ETC. 


jk d(m—n)(cos 6, cos Ot, —Sin 6, sin at, } 
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where 1 2 —-sin a, (A9) 
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1 J 
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-V. | cos -2—————sin$j)e — ?* f41-«^ dt 
zi bx JU 
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Let us now assume that b, >> 410 — be >> 10,50 that / b, 7 E = bx, then 


E 6 1 k.d(m- - 
Se on ev een) =| ap cost = T DK Tsin Oy )dt 
a k 
A 
_ 234) (kK d(m—n)sin 6, ——) 
= V2 e IK d(m-n)cos 9, SE 
(odtm nsn 2 
Ry 
- y AQ. non)eos ye _ _(kodsin dy. Rx) (A10) 
2J; e 
on k 
where f,,_, (k,dsin6, ‚re ,Rı) = — All) 
m-n (Ko Meth Rk) k,d(m—n)sin 0, Tp 
Rx 


and J, =Bessel’s function of the 1“ kind, of order 1. 


Substituting (AQ) into (A4), 
<V Ve inn = Vo? { <ap?> CFM MMOS GF (k dsind,,15,Ro) 


Sea Se oe eee O 


Therefore, c VANS SUV ue cte EX Te (A12) 
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k,d(m—n)sin oy ry 
Rx 

k d(m—n)sin Oy 1, 

Ry 


2 


jk ,d(m - n)cos 6, 


where C m E (A13) 


Einaly aD ZEV VS =< Ne eZ. SN ENN 
= VAN Cane en, > EZ 


EZ E EZ A Va 


= Vy [U + Ze Go Co) Zot + Ze (EE Cy) Ze* (A14) 
we as a . 
where & = — is the input signal to thermal noise ratio (A15) 
VN 
2 2 
SS > 
C= e is the input interference to thermal noise ratio (A16) 
VN 


B. EXPRESSION FOR S 


From (16) S = <A," V> 


where A, (reference signal) takes the form ao Vo D (A17) 


then S =<a V, e KGst*B,) (Ze [a eX st B.) Ve Y aye UB y 1 +N(0) > 
] 


= Vo Ze < ao” (Vor) + E (a, ay e PF - P9 y, 9o) (A18) 


+ Vo < a, e X9» B Nt 


< y,” >, E eJKod(n—-1) cos $, > 
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= V, eJKodin-1)cos0,r _ (kdsindutz.R ) (A19) 


Lastly, <a, a, Pk P9 V.%> = <a, a, > <elPk- PY s VS =0 (A20) 
where < a, a. > = 0. 
Therefore, using (A20), (A23) and (A24), 
S  =Vo <a >Zel fur(Kod, ro, Ro) ] 


= x ee Lat. (A21) 


where 


— 2), (k,d(m - D)sin 6, = 
[C; Be Luz cas Goo ee Ole (A22) 


(k ,d(m-—1) sin Do I 
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APPENDIX B : EXPRESSIONS FOR P, AND P, 


From (32), we have 


P, = < (wH Z. A, V4) (wB Z, A, V4) - 

Loca V eV VE E 

= UN E, wP Z, C, Zw (Bl) 
From (33), we have 
Pr =< (we ZL. Ax Vi’) (we Z. Ak NS > 


: 
mc) c wc Ne ew 


= VE WE Z ECO ZE y (B2) 


where Vy = root mean square noise voltage 


Es , & as defined in (A16) and (A17) 
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